Cluster Luttinger liquids of Rydberg-dressed gases in optical lattices 
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We investigate the zero-temperature phases of hard-core bosonic and fermionic gases confined to 
one dimension and interacting via a class of finite-range soft-core potentials. Using a combination 
of analytical and numerical methods, we demonstrate the stabilization of critical quantum liquids 
with qualitatively new features with respect to the Tomonaga-Luttinger liquid paradigm. These 
features result from frustration and cluster formation in the corresponding classical ground state. 
Characteristic signatures of these liquids are accessible in state-of-the-art experimental setups with 
Rydberg-dressed ground state atoms trapped in optical lattices. 
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Frustration is responsible for a variety of phenomena 
in many-body physics, such as the emergence of gauge 
symmetries |Tj, macroscopic ground state degeneracies 
in strongly-correlated quantum spin and fermionic mod- 
els [5] and the glass phases 0. While these effects per- 
tain to two and three dimensions, in one-dimension (ID) 
frustration usually leads to gapped states in either spin 
or fermion systems [T], while gapless quantum systems 
seem to invariably fall within the universality class of 
Tomonaga-Luttinger (TL) liquids [IHH]- There, accord- 
ing to Haldane's construction 0, the low-energy behav- 
ior is described by a free bosonic theory [9j, where the 
asymptotic decay of correlation functions in the micro- 
scopic model can be predicted in terms of just a few 
effective parameters, such as the Luttinger parameter 
K = vnC, with v and C the sound velocity and com- 
pressibility of the liquid [6^. While this can be proven 
for certain integrable systems, a fundamental question is 
whether the TL paradigm persists for any kinds of local 
or non-local interactions. 

Non-local (dipolar, van-der-Waals) interactions can be 
now realized experimentally with ultracold gases of po- 
lar molecules [101 E], groundstate magnetic [H] and 
Rydberg-excited atoms [I3HI6] confined to optical lat- 
tices [171 in]- In particular, in two and three dimen- 
sions soft-core interactions in gases of groundstate atoms 
weakly dressed with a Rydberg state using laser light 
[TOl [20] have been shown theoretically to generate ex- 
otic phenomena, such as, e.g., supersolid states of matter 
[2TI - [27] . Quantum phases of ID soft-core systems have 
remained so far largely unexplored. 

Here, we show that cluster formation due to soft-core 
potentials [551 US] on a ID lattice provides a natural 
realization of frustrated systems with infinite ground- 
state degeneracy. This allows for the stabilization of 
quantum liquids with anomalous cluster-type correla- 
tions not captured by the TL model. We derive an ef- 
fective field theory describing such states by means of 
bosonization of effective degrees of freedom on an in- 
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Figure 1. Sketches of the classical ground states (i — 0) 
of Eq. ([l| in the liquid (a), crystalline (b) and cluster (c) 
phases, for re = 2, 3 and 4, respectively and particle density 
n — 1/4. Particles are drawn in the lattice as blue spheres 
(see text). Insets: numerical DMRG computations of density- 
density correlations g2{x) for V/t = 6 and L = 48 in the cor- 
responding phases. For re — 4 [panel (c)] particles form two 
kinds of blocks (or clusters) with Np = 2 and Np — 1 particles 
{A and B in the figure, respectively). In the quantum regime, 
this results in a cluster-type liquid with correlation functions 
different from a TL liquid [see Insets]. 



finitely degenerate ground state manifold and confirm 
the departure from the TL picture by means of numerical 
simulations based on the density-matrix-renormalization- 
group (DMRG) algorithm JUl IS]- Gharacteristic fea- 
tures of these anomalous cluster Luttinger (cL) liquids, 
are evident in correlation functions and collective exci- 
tations. cL liquids should be realizable in experiments 
for Rydberg-dressed ground state alkali atoms, where 
the presence of the optical lattice potential allows one 
to reach the strongly interacting regime. 

The relevant microscopic Hamiltonian for hard-core 



bosons or spinless fermions in ID geometry reads 

i i 1=1 

Here, h\{hi) are bosonic/fermionic creation (annihilation) 
operators at the site z, rii — b\bi, and t is the tunneling 
rate on a lattice of spacing a, which, from now on, we 
set as unit length of distances. In the following, p — n 
(p = 1 — n) for n < 1/2 {n > 1/2), with n the particle 
density. As shown below, the soft-core model of inter- 
actions in Eq. (IT]) emerges naturally in clouds of cold 
Rydberg-dressed atoms, where the strength V and range 
Tc of interactions can be independently tuned using exter- 
nal fields. In the following, we focus on the bosonic case, 
discussing the fermionic one at the end of the manuscript. 

Equation (IT]) has a rich phase diagram, resulting from 
the combination of the following features of soft-core in- 
teractions: (i) their Fourier transform has negative com- 
ponents, which in classical systems are associated with 
free-space cluster formation |55]; (ii) the competition of 
Tc with the critical length r* = l/p — 1 leads to frustra- 
tion for re > r*. These effects result in the formation of 
highly degenerate cluster-type ground states in the clas- 
sical limit (i = 0), as well as in novel phases in the quan- 
tum regime {t > 0). Note that all correlation functions 
are here particle-hole symmetric, such that the phases for 
filling n and 1 — n are identical. 

Fig. [T] shows sketches of the possible ground states for 
n = 1/4 in the classical limit t = 0. These phases are: 
a liquid (r^ < r*), a crystal {re = r*), and a cluster 
state {re > r*). The classical liquid [panel (a)] is given 
by the summation of all particle configurations with the 
constraint that inter-particle distances are larger than r^ 
while the crystal has a particle every re + 1 lattice sites 
[panel (b)]. The ground state in the cluster regime is in- 
stead highly degenerate. We find that its structure can 
be computed analytically using the following procedure: 
we recognize that (i) for re > r* particles and holes are 
grouped into blocks of type A and B consisting of Np 
and Np — 1 particles, respectively, followed by a number 
of re empty sites [Np = 2 in Fig. Iljc)]. Other block sizes 
are energetically unfavored. The total number of blocks 
M for a system of length L is M = L{1 — p)/re. (ii) 
Any exchange of two blocks does not change the energy 
of the system. The classical ground state then consists of 
all permutations of blocks A and B with the constraint 
that their absolute number is fixed (e.g [AABAAB], 
[ABAAAB], ...). In this case, the ground state degener- 
acy grows exponentially with d = Af !/[(A//3)!(2M/3)!]. 

The formation of the correlated block-structure above 
in the ground state is well captured by the static struc- 
ture factor S{k) = Y.i,, e^'''^"^'^52(^ - j)/L, with g2{l - 
i) — {ngrij) — n^ the density-density correlation func- 
tion. This is shown in Fig. l2|^a) for the case n = 3/4, 
Tc = 4 > r* and t — 0, where S{k) peaks at a A:-value 
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Figure 2. (a) Static structure factor S{k) at fixed r^ = 4 and 
several V/t for a chain of length L = 48. Black solid line; 
classical case {t = 0). (b) S{k) for V/t — 1.5 and re = 1, 
(black circles), 2 (red squares), 3 (blue diamonds), 4 (green 
up triangles) and 5 (yellow down triangles). 



ke < 7r/2 (solid black line). Numerical DMRG results in 
the same figure show that this peak persists for t ^ V, 
where the ground state is a strongly interacting quan- 
tum liquid. This signals a departure from the TL pic- 
ture, where fcc — t^ 1^ for any strength of interactions, as 
discussed below. Key questions are (i) what is the un- 
derlying mechanism and (ii) what is the effective theory 
connecting the classical and quantum regimes. 

Low-energy theory: For bosonic and fermionic gapless 
systems, the low-energy physics is usually well captured 
by the universality class of the TL liquid. In Haldane's 
formulation, the operators 6j are phenomenologically de- 
scribed by long-wavelength density and phase fluctua- 
tion fields, ^{x) and '&{x), satisfying \d{y) ,^ ip{x)\ — 
iTTS{x — y), and the density operator reads 



i{x) = [n — V(^(x)'/7r 



1 \ ^ 2ip{-Knx--:p' (x)) 



(2) 



with ip'{x) = nnx — ip{x)/2. The discreetness of par- 
ticle density is captured by the last factor in Eq. ([2|, 
which determines the behavior of correlation functions 
beyond the hydrodynamic approximation, as, for exam- 
ple, 52(2:) = 71^ + K/{2'Kx'^) -\- acos(27rna;)/x^^. Here, 
a is a non-universal constant and K is the Luttinger pa- 
rameter, which embodies the effect of interactions in the 
low-energy theory. The form of g2{x) implies peaks (if 
any) in the static structure factor S{k) at momenta com- 
mensurate with density excitations, independently of the 
form and strength of interactions. 

A quantitative example for a lattice case is shown in 
Fig. [^b), where DMRG resrdts for S{k) are plotted as 
a function of k for V/t = 1.5 and n — 3/4 and several 
values of r^ (a detailed discussion of numerical results is 
given below). Because of finite range interactions, a peak 
at fc = 7r/2 is formed for re <r* — 3, consistent with the 
TL picture. For re > r* , however, the peak shifts to 
lower values k < 7r/2, incompatible with a TL liquid. 

These observations require a modification of the ef- 
fective theory to incorporate inhomogeneous cluster-type 
density distributions. In the low-density limit ren ^ 1 



and n <C 1, this is achieved by first defining the density as 
nix) = Y.m=i fixm)S{x - Xm)- Here f{xm) € {1, 2} la- 
bels the cluster size (here, clusters with one and two par- 
ticles). We then introduce a collective cluster counting 
field ipci{x), so that the density operator can be rewrit- 
ten as n{x) = Vificiix) J2m=i f{xm)S[ipci{x) ~ nm]. The 
field ^ci{x) differs sharply from its counterpart ip{x) in 
the standard TL theory, in that it takes integer values 
at the position of clusters and not of individual particles, 
as exemplified in Fig.^b). This opens the way to novel 
collective excitations, as shown below. 

After introducing ip'^i{x) — {2Trnax — ipci{x)l'2), with 
(J = M/Np the cluster density, the corresponding func- 
tional dependence of n(x) on (p'^i{x) reads 



n{x) = n WipciixY 



E 



akC 



2ik[<^ci (x)' —nnax] 



, (3) 



with Ofe constants. Similar to TL liquids, in the crit- 
ical regime the effective low-energy theory is here ex- 
pected to be a c = 1 conformal field theory (correspond- 
ing to a compactified boson model); however, Eq. ^ 
implies that physical observables will differ from those 
of TL liquids. As an example, g2{x) now reads 32(2;) — 
n^ -t-ai/x^ -I- [a2 cos(2TTnax)]/ x^ , with ai, a-i and 7 con- 
stants. Here, unlike TL liquids, the relation a\ = 7/(47r) 
is not satisfied, preventing a direct estimate of the de- 
cay of correlation functions from spectral properties by 
means of, e.g., level spectroscopy techniques [32] ■ 

The position fee of peaks in S'(fc) now depends on the 
interplay between density and interactions, and is pre- 
dicted to be fcc — 2mTa = 27r(l — p)/rc, with kc < n/2 
(see also above). As an example, for r,, = 4: and p — 1/4, 
kc = 37r/8 in perfect agreement with classical results 
[see Fig. [2ja)]. The small deviation of the S = 2-k/L of 
the theory from the DMRG results will be subject of fu- 
ture studies. Appreciable differences occur for correlation 
functions and spectral properties. We note that the pro- 
cedure above allows one to retain high-momentum fea- 
tures due to interactions in the low-energy theory (here 
corresponding to cluster formation) , and can be in princi- 
ple applied to other problems where the classical ground 
state presents infinite degeneracies. 

Numerical results: We investigate the phase diagram 
of Eq. (IT]) by performing numerical DMRG simulations 
for a given particle density n as a function of V/t and re- 
in order to avoid boundary effects, which can be signif- 
icant due to the finite range of interactions, we perform 
finite-size scaling in systems as large as L = 64 sites with 
periodic boundary conditions, keeping up to 1800 states 
per block and performing 5 sweeps. This ensures numer- 
ical errors of order 10~^ or less. 

An example for n = 3/4 is shown in Fig.pfa). There, in 
the considered range of parameters, for sufficiently small 
interactions, the low-energy phase is shown to be always 
a TL liquid, independent of r^. This is clearly seen for. 
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Figure 3. Left panel: numerical zero-temperature phase dia- 
gram of Eq. ([l| for density n = 3/4. The following phases 
are indicated: Luttinger liquid (red circles), crystal (blue 
squares), cluster Luttinger liquid (yellow circles). Green tri- 
angles define a phase which shows an anomalous entangle- 
ment entropy growth (see main text). Right panel: Configu- 
rations allowed for the quantum field (/'(a;): in the TLL sce- 
nario (top), v'(a;) takes integer values at the position of each 
particle (red circles). In the cluster Luttinger liquid (bottom), 
ifici describes the cluster density, and not the particle density. 
This points to a different nature of collective excitations. 



e.g., Tc = 2 < r*, where the system is equivalent to a XXZ 
model with next-nearest-neighbor interactions and finite 
magnetization. For re = r* = 3, instead, the TL liquid 
turns into an insulating Mott phase with a hole every 
four sites [blue squares, see also sketch in Fig. fl|b)] via a 
quantum phase transition of the Berezinskii-Kosterlitz- 
Thouless (BKT) type around V/t ~ 4, as discussed 
in Ref. 132 1. The most interesting behavior occurs for 
fc/i"* > 1 and V/t > 1, where the TL turns into a cL 
liquid. Below we analyze in detail this regime, and com- 
ment on the region with V/t ^ 1 and large r^/r* > 5 
(green triangles) towards the end of the paper. 

We characterize the phase diagram by focussing on: i) 
S{k); ii) the asymptotic decay of (72 and of the one-body 
density matrix Bij — {h\hj)\ Hi) the excitation spectrum 
Q.{k) = E}./ S{k) in the Feynman approximation, with 
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(6l6.+i + /j.c.)|0). (4) 



VL{k) is an upper bound to the exact spectrum |33j . 

According to the discussion above, differences between 
TL and cL liquids are evident in S{k) and g2- in partic- 
ular, in the TL phase S{k) exhibits a peak at kc — 7r/2, 
commensurate with the particle density, while peaks 
commensurate with the cluster structure emerge in the 
cL phase, the exact position depending on re- For suf- 
ficiently strong interactions V/t > 1, Figs. l2](a) and (b) 
demonstrate the crossover from TL to cL liquid increas- 
ing strength of interactions at fixed r^ > r* and as a 
function of Tc at fixed V/t, respectively. 

Correlation functions such as Bij decay algebraically 



both in the TL and cL phases, vahdating the assump- 
tion that both underlying field theories are c = 1 con- 
formal. However, as illustrated in Fig. l4Fa), a strong 
oscillatory behavior for B^j in the cL phase leads to the 
emergence of new peaks in the momentum distribution 
n{k) [Fig. [4Fb)]. This effect, which may be observed ex- 
perimentally, is predicted by the effective theory above 
and is reminiscent of the deformation of the Base surface 
observed in frustrated, two-dimensional bosonic models 
[Mj . Finally, as mentioned above, we note that estab- 
lished methods to extract K fail in the cL phase by giving 
non-consistent estimates in the re > r* region [32]. This 
further strengthen the departure from the TL scenario. 

The excitation spectrum is analyzed in Figs. Etc) and 
(d), where r2(fc) is plotted as a function of k for sev- 
eral values of r^ with fixed V/t = 3, and for increasing 
V/t with fixed re — 4, respectively. Panel (a) shows 
that finite-range interactions with re < r* induce a ro- 
ton excitation at A; = 7r/2, commensurate with density 
excitations, as expected. This is a precursor of the BKT 
crystallization at r^, — r* and V/t > 4. For r^ > r*, 
however, a new roton-type excitation appears at smaller 
k. This corresponds to the novel cluster-type excitation. 
This effect is evident in Fig. |4[d), where fl{k) develops 
cluster-type features at fc > 7r/4 already for weak inter- 
actions V/t = 1 (red circles), reflecting the structure of 
the classical ground state. 

Finally, in the green region with V/t 3> 1 and re/r* > 1 
in Fig. IsFa), we observe both (i) correlation functions 
consistent with a cL liquid and (ii) a significant increase 
of the entanglement entropy. This may point towards a 
critical phase with central charge c > 1. However, strong 
finite size effects prevent here any conclusive statements 
on this regime, which we plan to investigate in a future 
work. The phase diagram of Fig. ^a) holds for fermions as 
well, where the TL phase displays dominant density- wave 
order; hints on possible physics beyond TL for similar 
Fermi systems have been reported in Ref. ^35j . 

Experimental setup and ohservahles. The phases 
above may be realized using ultracold Rydberg atoms 
trapped in ID in an optical lattice. Here we consider 
^^Rb atoms where each ground state \g) is off-resonantly 
coupled to an excited Rydberg state |r) via, e.g., a two- 
photon laser with effective Rabi frequency $7 and red 
detuning A, with |A| ^ 17, [see Inset of Fig. [sla)]. 
We choose \r) = |43S'i/2), rj/(27r) = 1 MHz and 
A/(27r) = 10 MHz (we set the Planck' s constant h = 2n). 
The microscopic two-body interaction potential reads 

K-yd(i,j) = 1^ryd{l + [(« " J>Ac]^}-\ with V^ryd/(2^) - 

r2''/(167rA'^) ~ 125 Hz. The on-site hard-core constraint 
is enforced using, e.g., Feshbach resonances. The re- 
sulting cut-off radius re = (C6/2A)^/^ ~ 2.2/im is con- 
siderably larger than usual spacings a ~ 0.4 /im. We 
choose a lattice depth yo/(27r) > l3ER/{2Tr) = 50 kHz, 
such that the longitudinal trapping frequency 1^11/(277) = 
Er\/Vq/Er ~ 4.4 kHz > max{i, T^yd} limits the dynam- 
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Figure 4. (a) One body density matrix for r^ = 4 with 
interaction strengths V/t = 1,3,5. (b) Momentum distri- 
bution for the same parameters in (a), (c) Dispersion rela- 
tion at V/t = 3: for re = 2 (re = 5) a density roton mini- 
mum (cluster-like roton minimum) appear at the correspond- 
ing peaks in S{k) [see text and Fig. [2|b)]. (d) Dispersion 
relation for the same parameters in (a) and (b). 
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Figure 5. (a) Inset: Laser dressing scheme (see text) 
and resulting Born-Oppenheimer potential Vryd(i, j) for two 
Rydberg-dressed atoms (thick red line) compared to the 
square potential V of Hamiltonian (II]) (thin blue line) . Panel 
(a): g2{x) for several strengths Vryd/t and Tc = 4 [symbols as 
in (b)]. (b) S{k) for several Vryd/t in comparison to V/t = 6 
(black circles). 



ics to the lowest lattice band. Here, Ej^ — 7r^/(2ma^) is 
the recoil energy of atoms with mass m, while i/(27r) < 
45 Hz is the hopping rate in the lattice. Transversal trap- 
ping frequencies ujj_ ^ clI|| confine the system to a ID 
geometry. The small effective decay rate of the dressed 
ground state 7off — (J7/2A)^7r ~ I — 10 Hz ensures dom- 
inant Hamiltonian dynamics, with 7^ the bare decay rate 
of \r). Figure [5] shows that 52 (a^) and S{k) using Vrydii,j) 
are essentially indistinguishable from those of the poten- 
tial V of Eq.(fl]) for Vryd/t ^ I. This opens the way to 
the observation of the phases above in experiments. 

In conclusion, we have shown that soft-core potentials 
support quantum liquid phases beyond the Tomonaga- 
Luttinger scenario. These results, which derive from the 



existence of a macroscopically degenerate ground state 
in the classical limit, are of immediate interest for ex- 
periments with ultracold Rydberg atoms. Interesting ex- 
tensions of this work include the search for exotic phases 
in 2D, such as, e.g., frustration-induced Bose metals |3S] 
and emergent gauge theories [T]. 
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